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Abstract 

It is by now well established that, by means of the integration by part identities 
|l]] , all the integrals occurring in the evaluation of a Feynman graph of given topology 
can be expressed in terms of a few independent master integrals. It is shown in this 
paper that the integration by part identities can be further used for obtaining a linear 
system of first order differential equations for the master integrals themselves. The 
equations can then be used for the numerical evaluation of the amplitudes as well 
as for investigating their analytic properties, such as the asymptotic and threshold 
behaviours and the corresponding expansions (and for analytic integration purposes, 
when possible). 

The new method is illustrated through its somewhat detailed application to the 
case of the one loop self-mass amplitude, by explicitly working out expansions and 
quadrature formulas, both in arbitrary continuous dimension n and in the n — > 4 
limit. 

It is then shortly discussed which features of the new method are expected to 
work in the more general case of multi-point, multi-loop amplitudes. 



PACS ll.10.-z Field theory 

PACS ll.10.Kk Field theories in dimensions other than four 
PACS 11.15.Bt General properties of perturbation theory 



1 Introduction 



The usefulness of the integration by parts identities |I| for reducing dramatically the 
number of independent integrals to be evaluated in the calculation of multi-loop Feynman 
graph amplitudes, expressing them in terms of a small set of independent integrals (often 
called master integrals), is by now well established (see for instance Ref. 0). It is the 
purpose of this paper to point out that the same technique can be used for obtaining 
a linear system of first order differential equations in any of the external Mandelstam 
variables for the master integrals themselves 0. 

The new differential equations provide with a fairly complete information on the ana- 
lytic properties of the integrals. That information can be conveniently used, for instance, 
for the direct numerical evaluation of the integrals themselves, or for the analytic expan- 
sions around particular values of the external scalar variables or of the internal masses 
- and that by merely algebraic methods, without any attempt at carrying out the loop 
integrations occurring in the definition of the Feynman graph. The equations can also 
be exploited for attempting their explicit analytic integration (an achievement which is 
however expected to be possible only in a limited number of particular cases). As required 
by the integration by part algorithm, all the loop integrals are defined in n continuous 
dimensions; the n — > 4 limit is also easily worked out within the new approach. 

To give a flavour of the features of the method, the almost elementary cases of 1- 
loop vacuum (0-point) and self-mass (2-point) amplitudes is discussed in some details (in 
those cases the system of equations reduces actually to a single equation). While it is 
claimed that the approach can be extended, without significant loss of effectiveness, to 
more general multi-loop and multi-point amplitudes, it is obvious that the 1-loop 2-point 
case is too simple to be really significant, so that from this point of view this paper can 
be considered just an introduction to further work (which is in progress). In any case, it 
is hoped that the particular application of the new method presented in this paper might 
be of interest on its own, at least pedagogically. 

The paper is articulated in seven Sections. This is the first Section; Section 2 sets the 
notation and deals with the 1-loop vacuum amplitude; Section 3 establishes the differential 
equation for the 1-loop self-mass amplitude, which is the main result of the paper; Section 

4 discusses a number of expansions which can be worked out from the equation; Section 

5 recasts the equation as a quadrature formula; Section 6 is devoted to the n — > 4 limit; 
Section 7, finally, contains an outlook to the possible extensions of the method to multi- 
point, multi-loop amplitudes. 
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2 The 1-loop vacuum amplitude. 

We will use a Euclidean n-dimensional loop variable k, with volume element defined as 

d n k 



(27r) n - 2 ' 

the Minkowski loop variables are recovered by the Wick rotation 

d n k = I d n ~ l k dk n = -i( [ d n ~ l k dkn) =-i{ I d n k 



(1) 



o 

Mink \J / Mink 



The simplest 1-loop 0-point graph is shown in Fig.l; the associated amplitudes are 

, p \ f d k 1 / \ 

T[n,m ,—a) = / - — r — - — — ; r— , (2) 

where a is a positive integer. 




Fig.l: The 1-loop vacuum graph. 
An alternative dimensionless definition is 



r d n k 1 

T { n, {m /^-a)=^l 1 ^ WT ^, (3) 

where /i is a scale with the dimension of a mass. For simplicity, we will however use in the 
following Eq.(0), corresponding to Eq.(|3|) at fx = 1. The integration by parts identities 
are obtained starting from 

d n k d ( 1 \ _ 

(2tt)"- 2 dk„ V M (k 2 + m 2 ) a J ~ ° ' 

(the integral of the divergence vanishes for small enough n); carrying out the algebra, one 
obtains 

T{n, m 2 , -a) = ( 1 - - U ) -^T(n, m 2 , -a + 1) , (4) 



2 



or the equivalent equation 

T(n, m 2 , —a) = m 2 T(n, m 2 , —a — 1) . (5) 

With Eq.(H) all the integrals with a > 1 are expressed in terms of the "master integral" 
corresponding to a = 1, 

T(n, m 2 ) = T(n, m 2 , -1) = / > ( 6 ) 



while the continuation of Eq.(|5|) to a = gives 

(2tt 



V(/'./'r'.0)= / ttS^=0 



(the same applies to negative values of a ). 

T(n,m 2 ) is divergent in the n — > 4 limit; by using twice Eq.@ one gets 

T(n, m 2 ) = (n _ 2 ^_ 4) T(n, m 2 , -3) , (7) 

which provides with a convenient analytic continuation in n, as T(n, m 2 , —3) is finite at 
n = 4. According to the definition Eq.(^), scaling k by m and then expressing the loop 
momentum in n-spherical coordinates, <i n fc = m n K n ~ l dKdQ(n) one easily obtains 

2 oN m"- 6 fl(ii) 7 K^dK 



T(n,m 2 ,-3) 



(27r) n - 2 J (K 2 + lf 



where Q(n) is the n-dimensional solid angle. At n = 4, a simple explicit calculation gives 

T(4,m 2 ,-3) = -^; 

one can then write 

n—2 

T(n,m 2 ) = -^- -C(n), (9) 

(n — 2j(n — 4) 

where C(n) is a suitable dimensionless function of n , defined as 

C(n) = 8T(n, 1, -3) , (10) 
with the known limiting value at n = 4 

C(4) = 1 . (11) 
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In most applications C(n) can be kept as an overall factor; its expansion in n around 
n = 4, which can be obtained from Eq.s(|lCi|,|8|), is actually not needed if C{n) multiplies a 
finite expression in which the l/(n — 4) singularities have cancelled out. 

It is further to be noted that for n > 2 Eq.(^) implies 

T(n,0)=0. (12) 

For convenience of later use, let us also look at the 1-loop, 2-denominator vacuum 
graph of Fig. 2. 




Fig. 2: The 1-loop, 2kIenominator vacuum graph. 



Considering for simplicity only the case in which both denominators are raised to the first 
power, by partial fractioning one finds 

/ (2^ (k* + ml)(k» + mZ) = T^hrt ( T(n ' ml) ~ T(n ' ml) ) ' (13) 

showing that the amplitudes associated to the graph of Fig. 2 are just algebraic combina- 
tions of the amplitudes T(n,m 2 ,a) associated to Fig.l. 



3 The 1-loop self-mass amplitude. 

The 1-loop 2-point (self-mass) graph is shown in Fig. 3. 




Fig. 3: The p-4k>op self-mass graph. 
The associated scalar integrals are 
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af 2 2 2\ f ^ / \ 

b(n,m^,m 2 , — an, — a 2 ,p j = / 77-7 — 7 ttt 77 — 77 ttt 77 — , (14) 

v ' l! 2 r ; 7 (2vr)™- 2 (A; 2 + mf) Q i((p- A;) 2 + m^) Q 2 v ; 

where 011,0:2 are positive integers and the external vector p is also Euclidean ( p 2 > 
if p is spacelike; the timelike region is to be recovered by means of the usual analytic 
continuation). The integration by parts identities are derived from the equations 



d n k d ( 1 



(2tt)"- 2 ^ V M (k 2 + m\)^((p- k) 2 + m 2 2 ) a \ 



, (15) 



where the vector can be either p^ or k^, so that there are 2 identities for each choice 
of the pair ai,a 2 . Writing down explicitly the identities is elementary; as a consequence 
of the identities, it turns out that all the integrals defined in Eq . (|14"D with ai\ or a 2 larger 
than 1 can be expressed in terms of a single master integral corresponding to a.\ — <x 2 — 1, 
namely 

S^mlmlp 2 ) =jj^ 2 {k2 + ml){{ l_ ky + ml) ■ (16) 

Indeed, by writing explicitly the two identities corresponding to aci = a 2 = 1 in Eq.(|T5p, 
after some elementary algebra one finds 



S{n,ml,ml, -1, -2,p 



2\ 



[p 2 + (mi + m 2 ) 2 ][p 2 - (mi + m 2 ) 2 ] 

(n — 3)(m 2 — m\ — p 2 )S(n, m 2 , m\,p 2 ) + (n — 2)T(n, m 2 ) 



(17) 



The similar formula for S(n, mf, m 2 ., —2, —1, p 2 ) is then immediately obtained by exchang- 
ing nil with m 2 in Eq. (]T7|) ; by further using the relation 

S(n,m 2 i,m 2 2 , -(a x + 1), -a 2 ,p 2 ) = — -77^2 S ( n > m i> m 2> _Q, i> ~ a 2,p 2 ) , (18) 

at\ omf 

which follows at once from the definition Eq. (WM , one obtains the corresponding equations 
for higher values of ai,a 2 - 

Summing up: all the integrals S(n, m\, m 2 ., — ai, — a 2 ,p 2 ), Eq. (|l4|) , can be expressed 
as a combination of the self-mass master integral S(n, m 2 , m|,p 2 ) and the vacuum master 
integral T(n,m 2 ), defined in Eq . (|T6|) and Eq.(U), times the ratio of suitable polynomials 
in p 2 , m\ and m 2 — a well established and already known result. 
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We consider now the p 2 -derivative of S(n,m 2 ,m 2 ,p 2 ). We can obviously write 

2 & ml 2 2 2\ 1 & ml 2 2 2\ 



(19) 



By using the definition of S(n, m\ , m 2 ., p 2 ), Eq. (|TE|) , and then carrying out explicitly the 
p^-derivatives, one obtains 



r5 

at 2 2 2\ o 

Pn-^-S{n,m 1 ,m 2 ,p ) = -2 



1 



p 2 — p ■ k 



with the obvious identity 



p ■ k — — (k 2 + m 2 ) — ((p — A;) 2 + m 2 ) + p 2 — m\ + m\ 
the previous equation becomes 

p ll ^—S(n,m 2 1 ,m 2 2 ,p 2 ) = - (p 2 + m\ - m 2 2 )S(n,m\,m 2 2 ,~\,-2,p 2 ) 

- S(n,mj,m 2 2 ,p 2 ) + T(n } m 2 2} -2) ; 
by using Eq.(|j) and Eq.(|l7]), after minor rearrangements Eq.(|i]) can be written as 



d_ 

dp" 1 



Ot 2 2 2\ 

S^m^m^p 
+ (n - 3) 
- (n - 2) 



n-2 



2p 2 

p 2 + m 2 + wij 



rv 2 2 2\ 

S(n,m 1 ,m 2 ,p J 



S(n m 2 77i 2 p 2 ) 
[p 2 + (mi + m 2 ) 2 ][p 2 + (mi — m 2 ) 2 } ' ' ' 

(p 2 + m\ — m 2 )T(n, m\) + (p 2 — m\ + m 2 )T(n, m 2 ) 

2p 2 [p 2 + (mi + m 2 ) 2 ][p 2 + (mi — m2) 2 ] 



(20) 



That same equation can be derived, alternatively, starting from the definition of 
S(n,m 2 ,m 2 ,p 2 ) Eq.(|T6D and replacing the loop variable k by k/X, obtaining 

S(n,ml,m 2 2J p 2 ) = A 4 " n S(n, A 2 m 2 , A 2 m 2 , A 2 p 2 ) ; 

by acting on both sides with A 2 ^ one gets the familiar scaling equation 



d 







P Qp2 1 dm 2 



d 



mr 



dm 2 



+ 



4 — n 



at 2 2 2\ 

S(n,m v m 2 ,p 



which according to Eq.fllq) can also be written as 

d „, o os n — 4 

dp 2 



2 ^ ml 2 2 2^ 

p —S[n,m 1 ,m 2 ,p 



- ot 2 2 2^ 
-S(n,mi,m 2 ,p 



+ S(n, m x , m 2 , —2, — l,p ) + S^n, m x , m 2 , —1, — 2,p 
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by substituting Eq . fll7D and the similar formula for S(n, mf, mj, — 2, — l,p 2 ), Eq. (|20|) is 
recovered. 

Eq.(|20|) is the central result of this paper. It is an inhomogeneous, first-order differ- 
ential equation in the external variable p 2 for the master integral S(n,m 2 ,m 2 ,,p 2 ); the 
coefficients are ratios of polynomials in p 2 and the masses, the inhomogeneous term is 
known, consisting of the simpler vacuum amplitudes T(n, m 2 ) discussed at length in the 
previous section. As Eq . (p0|) is a first order equation in p 2 , it determines S(n ) m\ 1 m\ 1 p 2 ) 
once a single "initial value" for some value of the variable p 2 is provided; given that initial 
value, the numerical integration of the equation up to any desired value of p 2 can be easily 
carried out along any path in the complex p 2 -plane — and that for arbitrary continuous 
dimension n, the n — > 4 limit being just a particular case (see below). 

Eq.(|20|) can also be used for studying the analytic properties of S(n,ml,m 2 ,p 2 ), 
its value at specific values of p 2 or, more generally, its expansion around those values. 
Eq.flZOp is particularly well suited for the discussion of the properties of S(n,m 2 ,m 2 ,p 2 ) 
at the potentially singular points of the equation, namely p 2 = 0, p 2 = —(mi — m 2 ) 2 , 
p 2 = —(mi + m 2 ) 2 and p 2 — > oo; Eq.(^0]), in addition, can also be recast in the form of 
a quadrature formula which can then be used for the actual analytical integration in the 
n — > 4 limit. 



4 Expansions of the 1-loop self-mass amplitude. 

Let us start to discuss the point p 2 = 0; to that aim, Eq.(^Dj) can be rewritten in the fully 
equivalent way 

-^S(n,ml,m 2 2 ,p 2 ) = - ^— - ( S(n, mj, m 2 2 ,p 2 ) + T ^ m ^ — T (^ m 2) \ 

Op Zp y TTZ^ — 7772 / 

, n-3 ( 1 1 \ 2 2 2n 

2 \p 2 + (mi + m 2 y p 2 + (mi — m 2 ) z J 

n — 2 / mi — m 2 mi + m 2 

4mi(mi — ml) \p 2 + (mi + m 2 ) 2 p 2 + (mi — m 2 ) 2 

n — 2 / mi — m 2 mi + m 2 

4m 2 (mf — m 2 ) \p 2 + (m% + m 2 ) 2 p 2 + (mi — m 2 ) 2 

By inspection of the definition Eq.(|l^) (but without any explicit integration on the 
loop momentum k) S(n,m 2 ,m 2 ,p 2 ) is seen to be regular at p 2 = 0, i.e. the function and 




T(n,m 2 2 ) 



(21) 
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its derivatives are finite at that point. The apparent singularity of Eq . (|2l|) at p 2 = must 
therefore disappear: that amounts to the condition 



S(n, ml, m\, 0) 



T(n, ml) — T(n, m\) 



m\ 



(22) 



where T(n,m 2 ) is known, as already remarked. In other words, the equation itself gives 
the explicit value of S(n, m\, m 2 , 0), once the information of the regularity at p 2 = is 
provided. That same result, Eq.(|22]) may be also obtained by noting that at p^ = 
(implying p 2 = ) Eq . (|T6|) goes into Eq. (pT3|). 

Knowing that S(n, m\, m\,p 2 ) is regular at p 2 = 0, we can also expand S(n, m\, m\, p 2 ) 
in p 2 around that point 



S{n,m\,m\,p 2 ) = ^ S k (n, m\, m 2 2 )(p 



2\k 



(23) 



k=0 



where So{n,m\,m\) = S(n, m\, m\, 0) is of course given by Eq . (f2~!|) . By inserting the 
expansion Eq.(|23"D into Eq. (|21|) the coefficients Sk{n,m\,m%) are recursively determined; 
the first explicit values are 



Si(n, ml, ml) 



S 2 {n,ml,m 2 2 ) 



n{m\ — m?,) 3 



n — A)(ml + ml) [T(n, ml) — T(n, m^)) 
+2(n — 2) (mlT(n, ml) — m 2 {T{n, m^j 



n(n + 2)(m 2 — m 2 ) 5 
(n — 4) (^6(ml + ml) 2 — n(ml — ml) 2 ^j [T(n, ml) — T(n, m^)) 
-12(n - 2) (ml + ml) ( m 2 2 T[n, ml) — m^T(n, ml) ) . (24) 

Note that the knowledge of the coefficients of the expansion around p 2 = allows also the 
accurate evaluation of the r.h.s of Eq . (|2l|) for small values of p 2 when trying its numerical 
integration starting from p 2 = 0. 



A similar approach can be in principle worked out also at the point p 2 



-im 1 -m 2 



It is known, from the conventional study of the analytic properties of S(n,ml,m\,p 2 ) 
(which will not be repeated here), that S{n,ml,ml,p 2 ) is also regular at that particular 
value of p 2 (sometimes referred to as pseudo-threshold). Imposing the regularity at that 
point Eq.(|2"T|) gives 



S(n, ml, ml, —(mi — m 2 ) 2 ) 



n — 2 m 2 T(n, mf) — miT(n, ml) 



n 



2m\m 2 im\ — m 2 ) 



(25) 



It is a somewhat striking feature of this approach that the condition of regularity at 
p 2 = 0, determining completely S(n,m 2 ,mlp 2 ) through Eq. fl2T|) , which is a first order 
differential equation, implies also the regularity at p 2 = —(mi — m 2 ) 2 , a result which we 
were however unable to derive directly from Eq.(pl|) only, i.e. without explicit reference 
to the definition Eq.flTEp. 

For discussing the p 2 — > 00 limit, let us switch to the usual inverse variable y = l/p 2 , 
introducing the function 

Hn,ml,ml,y) = S (n,ml,rnl,~^ . (26) 

In terms of the new variable y, Eq.([H|) becomes 

d ti — 4 

— $(n, m 2 , m 2 , y) = — <£>(n, m 2 , m 2 , y) 

dy 2y 

n — 3 ( 1 1 \ , . 22\ 
+ — 1 + — 1 ^{n 1 m 1 ,m 2 ,y) 



(mi+m2) 2 y (mi—m2) 

n-2 ( m 1 -m 2 m x + m 2 \ mf 2 , 
+ ~a 7 — 5 9T 1 1 1 J- 

n-2 I m l -m 2 m 1 + m 2 \ 2 

+ 1 7 — 2 2T 1 1 T{n,m 2 ) . (27) 

4m 2 (mf - m 2 ,) + (mi+m2)2 y + (mi _ m2? ] 

The equation exhibits at y = the singular behaviour y"^"- -4 )/ 2 , so that the function 
must be expanded for small 7/ as the sum of two series corresponding to the singular and 
regular parts 

00 00 
$(n, ml ml y) = y'^ 2 £ ^ (n, m 2 , m\)y k + £ (n, m 2 , m 2 )// . (28) 

fc=0 fc=0 

According to Eq.s (|26| , p8|) , the asymptotic expansion of S(n, mf, m 2 ,p 2 ) for large p 2 is also 
given by 

00 1 00 1 

S(n, ml mlp 2 ) = (p 2 )^ 2 £ *«(„, m 2 , m 2 2 )-±- + £ $^(n, m 2 , m 2 )-^ . (29) 



By inserting the expansion Eq.(^) into Eq.(^) - or more directly the expansion Eq. fl29|) 
into Eq . (|2l|) , disregarding from now on the variable y - and imposing the equality of the 
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coefficients of the same powers of the expansion parameter one obtains 



$ { o\n, ml, mj) 
$1 (n, mf, m^) 

$ 2 J (n,m 1 ,m 2 ) 



0, 

T(n, ml) + T(n, ml) 



n -4 



m l ~~ m 2 ) ^( n ? m l) + 



n 



n 



ml -ml)T(n,ml) , (30) 



and 



<P\ (n, m 1 ,m 2 ) 

Q>2 m U m 2J 



(n — 3)(ml + m^^^^n, m 2 , m^) 



; (n - 3) (n - 4) (ml + m|) + 2(n - G^m^l ^\n,m\, m 2 2 ) (.31) 



Note that Eg . , by itself, is unable to fix completely all the coefficients of the expan- 
sion; the overall normalization of the singular part, <$>ffl(n, m\, m|) is in fact as yet not 
determined. But some more information on ^q\u, ml, ml) can be obtained from Eq.(|l7D, 
by writing it, according to Eq.fllB|), in the form 

d 



Of 2 2 2\ 

:S{n,m 1 ,m 2 ,p 



dm; 



[p 2 + (mi + m 2 ) 2 ][p 2 — (mi + m 2 ) 2 ] 
(n — 3) (rrii ~ m l ~ P 2 )S(n, ml, m 2 2 ,p 2 ) + (n — 2)T(n, ml) 



(n-2) 



p + m\ + m 2 ( 2 . 
2ml 



' rr( 2 

-1 (n, m 2 



(32) 



by inserting in it the large-p 2 expansion Eq. (p9|) and then equating the coefficients of the 
same powers of p 2 we obtain the additional equation 







dm 2 



^\n, ml, m 



; 



(33) 



as the same is true for the derivative with respect to m\, we can write as a conclusion 

^\n,ml,ml)=^\n) , (34) 

expressing the fact that the the still unknown leading coefficient $jj (n) is independent 
of the masses, but depends only on the dimension n. 

Due to the condition $o ( n > m\, ml) = 0, Eq. (|5(J|) and of the expansion Eq. fl2Tj|) , for 
n > 2 Eq.flU) implies 



lim (p 2 )- (n ~ 4)/2 S(n,ml,ml,p 2 ) 



p — >oo 



< S \n) , 



(35) 
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while for smaller values of n the dominant term of the expansion is <&[ (n, ml,m 2 ) /(p 2 ). 

The point p 2 = —(mi + m 2 ) 2 , finally, corresponds to the usual threshold and is there- 
fore another singular point for Eq.(|21"D and S(n, m\, m 2 ,p 2 ). It can be discussed, when 
required, along the lines sketched for p 2 — > oo. In analogy to Eq.(|29|) one can write the 
expansion around p 2 = — (mi + m 2 ) 2 as the sum of two series, 

oo 

S(n,m\,m\,p 2 ) = (p 2 + (m 1 +m 2 ) 2 ) M)/2 ^*[ s) (n,m|,m^)(p 2 + (m x + m 2 ) 2 ) k 

oo 

+ ( n ' m l m D(P 2 + K+ m2ffi (36) 



k=0 



by inserting the expansion in Eq.(^TJ), it turns out that the coefficients ^ (n, m 2 , m 
are completely fixed, while the ^^(n, m 2 , m 2 ,) are all proportional to a same constant, 
which can be taken to be ^q\u, ml, m 2 ) and is not determined by Eq.(|2T|) alone. The 
first values of the coefficients are 

lTr (r)/ 2 2n (n — 3)(n — 2) / , 2 2 \ 

*o n> m n m 2> = ~n 7 ; \ [m 2 T{n,m 1 ) + m 1 T{n,m 2 )) 

2mim 2 (mi +m 2 ) v ' 

^{n, m \,m\) = ( n ~ 2 -^-)^\n, m l, m 2 2 ). (37) 
\2(mi+m 2 ) z 8mim 2 J 

As already done for the large p 2 expansion, one can insert the expansion Eq.(|36[) and 
the explicit values Eq. (|37D in Eq . (j32|) , obtaining the equation 

-^-^\n,ml,m 2 2 ) = - -^-) ^\n,ml,m 2 2 ) . (38) 

om 2 \ 2m 2 m\ + m 2 J 

Its solution, symmetrized in m\,m 2 , is 

^\n,m\,ml) = { ™™ ] ^ 2 *P(n) , (39) 
[mi + m 2 ) n A 

where the still unknown coefficient Wq (n) depends only on n. 



2) 



5 The quadrature formula. 



The 1-loop self-mass case is simple enough to allow to recast the differential equation in 
the form of a quadrature formula. By putting 

(n-3)/2 



of 2 2 2 

S(n,m 1 ,m 2 ,p 



-(n-2)/2 



(P 2 

■o{n,m\,m\,p 2 ) 



p 2 + (mi + m 2 ) 2 ) (p 2 + 



m\ — m 2 : 



(40) 
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Eq.(^Cj) becomes 

9 , 2 2 2^ 

a(n,m 1 ,m 2 ,p ) 



1 



... L . , 2 (p 2 ) ( - 4)/2 [(p 2 + (m 1 + m 2 ) 2 ) (f + (mi - m 2 f 

(n — 2) (p 2 + m 2 — ml)T(n, m 2 ) + (p 2 — m 2 + m 2 )T(n, m 2 ) 
The above equation is equivalent to the quadrature formula 



-(n-l)/2 

(41) 



a(n,m 1 ,m 2 ,p) = a(n,m 1 ,m 2 ,p ) 

p 2 



n-2 



c/x x^ n 4 ^ 2 (x + (mi + m 2 ) 2 ) (x + (mi — m 2 ) 2 ) 



-(n-l)/2 



(x + m 2 — m 2 )T(n, m 2 ) + (x — m 2 + m 2 )T(n, m 2 ) 



(42) 



where p 2 is some suitable value of p 2 . As S(n,m1,m 2 ,p 2 ) is regular at p 2 = 0, for n > 2 
Eq. (fl0|) implies a(n,ml,m 2 , 0) = 0, so that p 2 = is a natural choice. 

For m 2 = 0, by using Eq. (|I~2"D , one has 



n-2. 



a(n,m 2 ,0,p 2 ) = -- — -T(m,m 2 ) / rfx x (n " 4)/2 (x + m 2 )- (n - 2) ; 

o 

which becomes, by using Eq.(^|) and scaling x by m 2 within the integrand, 

p 2 /m 2 



(43) 



cr(n, m 2 , 0,p 2 ) 



C(n) 



2(n - 4) 



dxx (n - 4)/2 (x + l)- (n - 2) . 



(44) 



By comparing Eq.(|3~5[), Eq.(||]) and Eq.fl44]) a closed expression for $f } (n) is obtained 



$S s) (n) 



C7(n) 



2(n -4) 



. oo 

^- fdx x("- 4 )/ 2 (x + l)-("- 2 ) 



(45) 



6 The n — > 4 limit. 

It is known that in the n — > 4 limit 1-loop amplitudes are singular at most as l/(n — 4); 
from Eq.(Q) one has indeed, for n — > 4, 



T(n, m 2 ) = m 2 



CM 1, 1 ^. 

- v ; - + - mm - - + Ofn - 4) 
2 (n - 4 2 4 V ; 



(46) 
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where C(n) is given by Eq.s( |IT^|Tl~|) . If mm, is perceived as disturbing, replace Eq.(H) with 
Eq.([3|) and lnm goes into In (m//i). 

Accordingly, S^n, m 2 , m 2 .,^ 2 ) can be expanded in the n — > 4 limit as 

S(n,m 2 ,m 2 ,p 2 ) = < ^^-S (m 2 l} m 2 2 ,p 2 ) + Si(m 2 , m 2 ,p 2 ) . (47) 
n — 4 

Eq. (p0|) is an identity in n; therefore, it implies a set of independent equations for the 
various terms of the expansion in (n — 4). The part singular in (n — 4), i.e. the coefficient 
of C(n)/(n — 4) , gives 

-S^m^m^p 



dp 2 ' ' yp 2 2(p 2 + (mi + m 2 ) 2 ) 2(p 2 + (mi — m 2 ) 2 ) / 

S (m 2 ,m 2 ,p 2 ) + £) . (4* 



On the other hand, we can expand S(n, m 2 , m|, 0), Eq. (|22"D ; that gives, in the notation of 

S (m 2 ,m 2 ,0) = -i (49) 

and 

2 2 1 m?lnmi-mpnm 2 , . 

Si(m 1; m 2 , 0) = +- — — = 2- . (50) 



(Again, In mi becomes ln(mi//x) etc. if the r.h.s. of Eq . (|i~4]) is multiplied by fi A n ). 
With the input Eq.(f49"D, Eq. (|48|) gives, for any p 2 , 

So(m 2 ,m 2 ,p 2 ) = -- , (51) 

independent of p 2 and of the masses. 

For discussing the finite part of Eq. fl2"0|) in the n — > 4 limit, it can be convenient to 
write Si(m 2 , m 2 ,p 2 ) as 

Si(m?,mf,p 2 ) = Si(ml,ml,0) + S(m\,m 2 2 ,p 2 ) , (52) 

where Si(m 2 , m 2 ., 0) is known from Eq.(|50|). Note that Eq . ([52"D , due to Eq . (pi]) , amounts 
to 

S(ml,ml,p 2 ) = lim (S(n, m 2 , m 2 ,p 2 ) - S(n, m 2 , m 2 , 0)) , (53) 
or, in equivalent way, 

S(n, m 2 , mj, p 2 ) — ^ S^n, m 2 , m 2 ,, 0) + S{m\ } m 2 ,, p 2 ) . (54) 
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From Eq.(pC|) - or, rather, from its equivalent form Eg . (pT|) - one then obtains for 
S(m\ , m 2 ,p 2 ) the eguation 

(9 1 

u ~o/ 2 2 2\ 1 Tr/ 2 2 2\ 

—S{m 1 ,m 2 ,p ) = S(m 1 ,m 2 ,p ) 

+ ^ 1 2 i ? 1 , \2 + 2 i 7 \2 I ^( m ?> m 2,P 2 ) 

2 \p z + {mi + m2) P + — va-i) J 

m 2 mi , /mA 1 



m 2 — m 2 \m 2 ) [p 2 + (nil + m 2 ) 2 ][p 2 + (nil — 1712) 2 } 

1/1 1 \ 
I (55) 

8 \p 2 + (mx + m 2 ) 2 p 2 + (nil — m 2 ) 2 ) 

Eg.(|55|) can be used, in the same way as Eg.(^), for obtaining the coefficients of the 
expansions of S(m 2 ,m 2 ,p 2 ) around particular values of p 2 or a guadrature formula. 

The expansion at p 2 = is straightforward. By writing 

S(mlm 2 2 ,p 2 ) = Si(mlm 2 2 )p 2 + 0((p 2 ) 2 ) , (56) 

Eg. (1)3) gives 



S, ( m 2 m 2 ) - m i m ' ln (rni\ _ m 2 + m 2 
bx{m ^ m2 > ~ 2{m\-mlf \mj 8(m 2 - m 2 2 ) 2 ' l&7j 

which is just the n — > 4 limit of S\[n, nil, m 2 ), Eg. (|24|) . 

More interesting is the n — > 4 limit of the p 2 —>■ 00 expansion. Let us start with 
Eg. (p9|) ; for n — > 4 one has 

(p2) (n-4)/2 ^ x + ^ ^ p 2 (5g) 

(see the comments to lnmi in Eg.(50) for the appearance of lnp 2 ). As it is known that 
one loop amplitudes have at most simple poles in (n — 4), the leading coefficient ^^(n), 
Eg.(^), can be expanded as 

$« (n) r^|^M 0o + 01) (59) 

n — 4 



where C(n) is again given by Eg.s(10],ll]) and 4>o, 0i are as yet unknown constants. In 
the n — > 4 limit the large p 2 expansion of Eg.fl29|) then reads 
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S(n,ml,m%,p 
+ 



2 2\ ri^A C(n) ,.,.1,1 2 

" — ► T0o + 0i + -^olnp 

n — 4 2 



2 2 JC(n) 20 o + l 
K+m 2 ) 2 

H — (m 2 In mi + m 2 In m2 



+ 0o + 0i ~\ + ln p2 
1 J 1 \ 



+ C 



(60) 



It is however known from Eg . s (|47|j5T[) that the pole in (n — 4) is independent of p 2 and of 
the masses; that implies 

0o = -\ , (61) 

so that Eg. (|60|) becomes 

C(n) 1 



of 2 2 2\ 

S{n,m 1 ,m 2 ,p ) — ► 



2(n-4) +01 -4 lnp 



+ 



(0i " ^ J K + ™ 2 ) - 4^1 In I ^2 I - tW 2 In 



P 

ml 



- + o J- 

2 Up 2 ) 2 , 



, (62) 



where the constant 0i is still unknown. 

Let us now look at the quadrature formula. We have already seen that the replacement 
Eq.([H]) transforms Eq.(pO|) into the quadrature formula Eg. fl4"2|) . One can similarly put 



S(ml,mj,p 2 ) = -J(p^ + K + m 2 ) 2 ) (p 2 + (m x - m 2 ) 2 ) • cf(m 2 , m 2 2 ,p 2 ) 

pz V 



(63) 



according to that definition, cr(m 2 ,m 2 ,0) = 0. Eq. (|55|) then becomes the quadrature 
formula 

r>2 



—/ 2 2 2\ 

a(m 1 ,m 2 ,p ) 



dx 



{(x + [mi + m 2 ) 2 ) (x + (m x — m 2 ) 2 )] 2 

1 , o 2 . m\m\ . /mi 

— -x[x + m 1 + m 2 ) H 5 ^ In — 

4 mf — m 2 \m2 



(64) 



An elementary integration gives the known j|] explicit analytic expression 
£(m 2 ,m 2 ,p 2 ) 



1 m? + m\ (m\ 

+ — — k rrr m 



4 4(m 2 — m 2 ) Vm 2 



+ 



4p2 



where 



{p 2 + (mi + m 2 ) 2 ) (p 2 + (mi — m2) 2 ) lntiQo 2 ) + (m 2 — m 2 ) ln 
wp 2 + (mi + m 2 ) 2 - -y/p 2 + (mi - m 2 ) 2 

w(p ) — 



mi 
m 2 



p 2 + (mi + m 2 ) 2 + wp 2 + (mi - m 2 ) 



,(65) 



(66) 
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to 



For real and positive p 2 , u(p 2 ) is also real and positive, ranging from u(+oo) 
u(0) = m 2 /m 1 (if m 2 < toi). For real and negative p 2 > -(mi — m 2 ) 2 , u(p 2 ) remains real 
and positive, with u(—(mi — m 2 ) 2 ) = 1; at p 2 = —(mi — m 2 ) 2 , further, Eq.fl6"5|) reads 



S(ml,m 2 , -(mi - to 2 ) 2 ) 



1 



m\m 2 



2(rrfi 



to| 



ln('^ 

m 2 



(67) 



o, . n ln ( — ) > ( 6J 

2(toi — m 2 ) \m 2 



From Eq. (|2q) and Eq.(|22|) a simple explicit calculation gives 
lim (S(n, m?, m^, — (mi — m 2 ) 2 ) — S(n, ml, ml, 0)) = - + ' 

so that Eq. (|53|) is satisfied at p 2 = —(mi — m 2 ) 2 . 

When p 2 varies from the pseudo-threshold —(mi — m 2 ) 2 to the threshold —(mi+m 2 ) 2 , 
one can put 

\jp 2 + (mi - m 2 ) 2 = i\J-(p 2 + (toi - m 2 ) 2 ) ; (69) 

ii(p 2 ) then varies on the complex unit circle, u(p 2 ) = e t ^ p2 \ with (f>( — (mi — m 2 ) 2 ) = 
and 4>( — (mi + m 2 ) 2 ) = — tt, so that in Eq.(|65D \nu(p 2 ) becomes i, the imaginary unit, 
times an arctangent, the square root in front of it is also imaginary and S(m 2 ,m 2 ,p 2 ) 
remains real (as S(m 2 ,m 2 ,p 2 ) is analytic at the pseudo-threshold, changing i into — % in 
Eq.(|69|) does not change its value). For p 2 < —(mi + m 2 ) 2 , finally, u(p 2 ) is negative, 
ranging from —1 at p 2 = —(mi + m 2 ) 2 to at p 2 = —oo; correspondingly, the square root 
of Eq.(|65|) is again real, while lnu(p 2 ) becomes complex, acquiring an imaginary part m 
whose sign however can be determined only by giving an infinitesimal imaginary part to 
p 2 (if Imp 2 < 0, the imaginary part is m\ 

Further, on account of the explicit analytic expression Eq.(^), in the p 2 — > oo limit 
Eq.(Bl) becomes 



of 2 2 2 

5(n,TO 1 ,TO 2 ,p 



n^4 



C(n) 1 1 2 

v ' r H Inp 2 

4 2 4 y 



m\ + m A 2 + m\ ln I ^—z I + m\ ln ( 



2(n 
p 2 ^ 



m 



m$ 



1 



(p 



2\2 



(70) 



the comparison with Eq.(|62|) gives 

<h = \- (71) 

Eq . s (|6TUTT|) can be obtained also by the explicit integration of Eq. (|4~5|) in the n — > 4 limit. 

Similarly, one can expand Eq. (|65|) around the physical threshold p 2 = —(mi + m 2 ) 2 , 
obtaining, in the notation of Eq . (p9|) , 



^(4) = — 



7T 



(72) 
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7 Outlook. 



As already said in the Introduction, the 1-loop self-mass is too simple for really showing 
the potential of the new method. It is nevertheless possible to try to assess which features 
of the method can be reasonably expected to apply also to less elementary cases. 

To start with, it is apparent from the derivation that a linear system of first order 
differential equations can be established, in full generality, for multi-point and multi-loop 
graphs as well, as the integration by part identities can be written for any number of 
external lines and loops. 

In the case of more general multi-point amplitudes, depending on a set of independent 
external vectors Pi ifJr , one has simply to generalize Eq. (|T9|) by considering all the various 
Pi ujt— combinations and then work out some kinematical algebra for extracting the 
derivative with respect to any given external scalar variables. A first application to the 
1-loop vertex is underway []5J. 

In the general multi-loop case several independent master integrals are expect to ap- 
pear, so that rather than a single first order equation there will be a linear system of first 
order differential equations, whose coefficients however will be in full generality ratios 
of suitable polynomials in the dimension n, the masses and the external scalar variables 
(work on the 2- loop "sunrise" self- mass graph is in progress ||). In the general case, the 
equations will involve, besides the master integrals for the considered original Feynman 
graph, a number of "inhomogeneous terms" corresponding to the master integrals for all 
the related, simpler Feynman graphs (with the same number of loops) obtained by re- 
moving a denominator from the original graph; those simpler graph amplitudes are to be 
considered known — or else can be studied recursively by the same algorithm. 

In the general multi-point, multi-loop case and for arbitrary values of the masses simple 
quadrature formulae and closed analytic results in terms of known functions cannot be 
expected to exist, but the equations will be anyhow a powerful tool for studying the value 
of the amplitudes for particular values of the occurring masses as well as for attempting, 
in particularly simple cases, the analytic integrations. 

The equations are particularly well suited for the expansion in (n— 4); as the equations 
are identities in n, the n — > 4 limit can be worked out at once for the whole multi-loop 
amplitudes, without any reference to the n — > 4 limit of the inserted loops, so that any 
overlapping divergencies problem is avoided. 

Once a system of equations for quantities finite in the n — > 4 limit has been obtained, 
its numerical integration is straightforward, for virtually any number of loops and external 
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lines, provided suitable starting points are given. It is reasonable to think (but also 
confirmed by some preliminary work, [H]) that the absence of the would-be kinematical 
singularities, such as those at the p 2 = for the self-mass, will provide with useful 
information for the evaluation of those starting points. 

The equations will also be very useful, in general, to provide any kind of required 
expansion of the graph amplitudes, for instance for asymptotically large values of the 
scalar variables or at the physical thresholds, in terms of very few constants; even if those 
constants are not expected to be fixed by the equations themselves, the information which 
they contain can be exploited to make easier the calculation of those missing values. 

Acknowledgments. The author is pleased to thank M. Caffo, H. Czyz and S. Laporta 
for several clarifying discussions. 

References 

[1] K.G. Chetyrkin and F.V. Tkachov, Nuc. Phys.B192, 159(1981); F.V. Tkachov, Phys. 
Lett.BlOO, 65 (1981). 

[2] S. Laporta and E. Remiddi, Acta Phys. Pol. B28, 959 (1997). 

[3] A differential equations method, based on mass derivatives, has been proposed by 
A.V. Kotikov, Phys. Lett. B254, 158 (1991). In that approach, amplitudes with a 
single non vanishing mass m are expressed as a suitable integral of the corresponding 
massless amplitudes, which are taken as known. 

[4] See for instance Eq.(B.l) of M. Bohm, H. Spiesberger and W. Hollik, Fortschr. Phys. 
34, 687 (1986). 

[5] M. Caffo and E. Remiddi, in progress. 

[6] M. Caffo, H. Czyz and E. Remiddi, in progress. 



18 



